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ABSTRACT 

I present a generalized power-law diagnostic that allows to identify the presence of active galactic 
nuclei (AGN) in infrared (IR) galaxies at z > 1, down to flux densities at which the extragalactic 
IR background is mostly resolved. I derive this diagnostic from the analysis of 174 galaxies with 
5^(24 jam) > 80 /iJy and spectroscopic redshifts z spec > 1 in the Chandra Deep Field South, for which 
I study the rest-frame UV/optical/near-IR spectral energy distributions (SEDs), after subtracting a 
hot-dust, power-law component with three possible spectral indices a = 1.3, 2.0 and 3.0. I obtain 
that 35% of these 24 /im sources are power-law composite galaxies (PLCGs), which I define as those 
galaxies for which the SED fitting with stellar templates, without any previous power-law subtraction, 
can be rejected with > 2<j confidence. Subtracting the power-law component from the PLCG SEDs 
produces stellar-mass correction factors < 1.5 in > 80% of cases. The PLCG incidence is especially 
high (47%) at 1.0 < z < 1.5. To unveil which PLCGs host AGN, I conduct a combined analysis of 
4Ms X-ray data, galaxy morphologies, and a greybody modelling of the hot dust. I find that: 1) 77% 
of all the X-ray AGN in my 24 /im sample at 1.0 < z < 1.5 are recognised by the PLCG criterion; 2) 
PLCGs with a = 1.3 or 2.0 have regular morphologies and Td us t ^ 1000 K, indicating nuclear activity. 
Instead, PLCGs with a = 3.0 are characterised by disturbed galaxy dynamics, and a hot interstellar 
medium can explain their dust temperatures Td us t ~ 700 — 800 K. Overall, my results indicate that 
the fraction of AGN among 24 jam sources is between ~30% and 52% at 1.0 < z < 1.5. 

Subject headings: infrared: galaxies - galaxies: evolution - galaxies: high-redshift 



1. INTRODUCTION 

The study of the individual dust-obscured sources 
that make the extragalactic infrar ed (IR) background 
(jPuget et al. 19961 : iDole et al. 20061 ) has made enormous 
progress over the last de cade, since the lau nch of the 
Spitzer Space Telescope (IWerner et al. 2004}). and suc- 
cessively the Akari Telescope ({Muraka mi et al. 2QQ7T) and 
the Herschel Space Observatory (jPilbratt 2003[ ). In par- 
ticular, the scientific output of these missions has re- 
vealed the importance of powerful, dust-obscured star- 
formation and nuclear activity in shaping galaxy evo- 
lution at high redshifts [z > 1). This activity was 
dominated by luminous and ul tra- luminous infrared 
galaxies (LIRGs and ULIRGs; iLe Floc'h et al. 20051 : 
iCaputi et al. 20071 ). which had a main role in the 
global star formation h istory of the Universe (e.g. 
iHopkins fc Beacom 20061), and the process of massive 
galaxy buildup at z > 1 (jCaputi et al. 2006aT ). 

The availability of multi-wavelength ancillary data 
from deep galaxy surveys has been crucial to investi- 
gate the presence and properties of IR galaxies at high 
z. Most redshift estimates and the derivation of other 
parameters, such as the galaxy stellar mass, rely on the 
fitting of spectral energy distribution (SED) templates 
to broad-band photometry that traces the galaxy rest 
UV/optical and near-IR light (A res t ^ 3 jam). The deriva- 
tion of reliable values for these galaxy parameters re- 
quires a proper wavelength coverage of the photometric 
data, and also considering the possible variations that 
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the galaxy SEDs may have, through the choice of suffi- 
ciently representative galaxy templates, and the analysis 
of departures of these SEDs from models of pure stellar 
evolution. 

A long-standing problem regarding the composition 
of the extragalactic IR background is understand- 
ing to which extent dust-obscured active galactic 
nuclei (AGN) are part of the IR galaxy popula- 
tion at high redshifts. A direct way to reveal nu- 
clear activity is searching for X-ray detections (e.g . 
Rigbv et al. 20041: iPolletta et al. 2001 iFiore et "aT~ 2008: 



Svmeonidis et al. 2011c iMatsuta et al. 20121 ). but dust- 



obscured AGN can remain undetected even in typically 
deep X-ray maps. Spectral line diagnostics at opti- 
cal and IR wavelengths are also useful to recognise 
AGN, but they are usually limited to relatively small 
samples of bright IR sources (e.g. [ Armus et al. 2007; 
[ Saiina et al. 20071: [Caputi et al. 20 08: De sai et al. 2008 : 
Nar dini et al. 20081: iHernan-Caballero et ~ 2009: 

iPetric et al. 20111 : [Hwang et al. 20121 ). 

An alternative method that has proven to be quite ef- 
ficient to reveal AGN in large samples of dust-obscured 
galaxies is the u se of mid-IR col our-colour diagrams 
iStern et al. 2005b . based on photom 



(|Lacv et al. 2C 

etry at 3.6 — 8.0 /im, t aken with the Sp itzer Infrared 
Array Camera (IRAC; iFazio et al. 2004f ). The seg- 
regated locus that some AGN occupy in these dia- 
grams is related to the fact that pure AGN SEDs 
are characterised by a power-law shape at rest-frame 
optical/near-IR wavelengths. In fact, the presence of 
a pure IRAC power-law SED has been proposed by 
some authors as a criterion to select AGN-dominated 
galaxies among Spitzer 24 /im and IRAC-selected sources 
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(jAlonso-Herrero et al. 2003 iDonlev et al. 2007D . But 
AGN selection through IRAC colour-colour diagrams has 
well-known li mitations: it is co mplete only for bright IR 
sources (e.g. lLacv et al. 2 007). and the resulting sam- 
ples are typically co ntaminated by star-forming galaxies 
(jDonlev et al. 20121 ). This is especially the case at high 
redshifts, as the galaxy stellar emission is shifted into the 
IRAC bands, producing similar IRAC colours to those of 
AGN at lower redshifts. 

In addition to the identification, there is a sec- 
ond problem, which consists in understanding what 
fraction of the galaxy light is due to the AGN com- 
ponent at different wavelengths. At mid- and far-IR 
wavelengths (3 < A res t < 1000 /im), disentangling the 
star- formation/ AGN contributions is required to prop- 
erly derive the on-going, obscured star formation rates 
in IR galaxies. This issue has been tackled in different 
studies of LIRGs and ULIRGs at high z, by analysing 
mid-IR spect roscopic data, or m id- /far-IR broad-band 
data (e.g. [ Fadda e t al. 201(1 iGruppioni et al. 2010: 
Bar thel et al. 20121: lHanami et al. 2012 : 

Mullan eTet al. 20121 : iPozzi et al. 20121 ). 

In the rest-frame optical/near-IR (0.3 < A rest < 3 /im), 
a power-law component associated with an AGN dis- 
torts the light of the underlying stellar populations in 
the host galaxy, and may affect the derived galaxy stel- 
lar mass. For galaxies that are identified with X-ray 
luminous AGN, the SED fitting with stellar templates, 
and derivation of ste llar masses are usually avoided (e.g. 
iCaputi et al. 2006bT ). or attempted only after subtract- 
ing empirical AGN templates (e.g. iMerloni et al. 20Tot 
iMainieri et al. 201 If ). In all other cases, however, the 
derivation of stellar masses is of common practice, even 
when there is the suspicion that the rest-frame near-IR 
light could be partly contaminated by an AGN compo- 
nent (for example, in the case of an IRAC-band excess 
in the galaxy SED). Quantifying the importance of this 
effect in IR-selected galaxies is then necessary to under- 
stand the reliability of the derived host galaxy properties. 

In this work I present the analysis of the rest-frame 
optical/near-IR SEDs of 174 24 /im-selected galaxies with 
S u {24fim) > 80/iJy and secure spectroscopic redshifts 
z spec > 1. My aim is to investigate whether, and 
when, a power-law component makes a relevant con- 
tribution to the galaxy SED, i.e. the subtraction of 
a power-law component produces a significant improve- 
ment of the galaxy SED fitting with stellar templates. 
Hereafter, I will refer to these galaxies as power-law- 
component galaxies (PLCGs). Note that this criterion is 
more relaxed than the IRAC power-law shape imposed 
by other authors to select p otential IR AGN candidates 
(jAlonso-Herrero et al. 20061 ). and therefore, the sample 
analysed here should also include galaxies for which the 
AGN emission is less dominant over the underlying host 
galaxy light. 

This paper is organised as follows. In Section $2j I 
give details on the sample selection and compilation of 
multi-wavelength photometry. In Section JJJ I present 
the results of the SED analysis, and quantify the ef- 
fect on the derived stellar masses. Later, in Section 
I investigate the nature of the PLCGs at z > 1. I 
make use of ultra-deep X-ray data to obtain an inde- 
pendent diagnostic of nuclear activity among the 24 /im 



sources, and understand in which cases the significant 
SED power-law component is related to the presence of 
an AGN. I also analyse the PLCG morphologies, and 
derive the greybody temperatures that are necessary to 
explain the SED power law. In Section $5j I investigate 
how the latest proposed IRAC colour-colour criteria deal 
with the selection of PLCGs and other X-ray-detected 
AGN. Finally, in Section $6j I summarise my findings 
and present some concluding remarks. All magnitudes 
and colours quot ed in this paper are total and refer to 
the AB system (jOke fc Gunn 1 983). I adopt a cosmology 
with H = 70kms- 1 Mpc- 1 , Sl M = 0.3 and ft a = 0.7. 
All stellar masses refer to a Salpeter (1955) initial mass 
function (IMF) over star masses of (0.1 — 100) M . 

2. SAMPLE SELECTION AND MULTI- WAVELENGTH 
DATASETS 

The Great Obser vatories Origins Dee p Survey 
(GOODS) programme (jGiavalisco et al. 2004f ) comprises 
a wide range of deep galaxy surveys conducted with main 
astronomical facilities, including the Spitzer and Hubble 
Space Telescopes, the Chandra X-ray Observatory, and 
the largest optical/near-IR ground-based telescopes. In 
particular, the Spitzer images for the GOODS fields have 
been collected as part of the GOODS Spitzer Legacy Pro- 
gramme (PI: M. Dickinson), and include deep data from 
both IRAC an d the Multiband I maging Photometer for 
Spitzer (MIPS; iRieke et al. 200l . 

The 174 24 /im-selected galaxies analysed here have 
been extracted from the publicly released 24 /mi cata- 
logue of the GOODS-South (GOODS-S) field, which has 
a flux density limit of S v (2Ajim) = 80 /iJy, down to 
which this catalogue is highly reliable and completfl 
Sources down to this flux density li mit make ~ 70% of the 
extragalactic 24 /im background (Pa povich et al. 20041 : 
iDole et al. 20061 ). 

To decide which 24 /im sources would be part of the 
analysis sample, I applied the following two criteria: 1) 
the source should have a secure (good-quality flag) spec- 
troscopic redshift (z spec ) from any of the multiple spec- 
troscopic galaxy surveys of the GOODS-S field (Le Fevre 
et al. 2004; Szokoly et al. 2004; Teplitz et al. 2007; 
Vanzella et al. 2008; Popesso et al. 2009; Fadda et 
al. 2010; Kurk et al. 2013); 2) the redshift should be 
z spec > 1. I searched for spectroscopic counterparts of 
the 24 jam sources within a 1.5 arcsec matching radius, 
and considered only one-to-one identifications. The final 
sample contains 174 24 /im sources with z spec > 1, with 
similar numbers of galaxies at 1.0 < z < 1.5 and z > 1.5 
(89 and 85 sources, respectively). 

Restricting the SED study only to galaxies with se- 
cure spectroscopic redshifts ensures that all the results 
and conclusions in this paper are free from the uncer- 
tainties that are usually introduced by photometric red- 
shift determinations at high redshifts. On the other 
hand, the multiple GOODS-S spectroscopic datasets 
correspond to a wide variety of selection criter ia, and 
even include some IR spectra (Fad da et al. 2010l h so the 
sub-sample of 24 /im galaxies with spectroscopic red- 
shifts is reasonably representative of the entire sample 
of 5' i/ (24/im) > 80/iJy galaxies at z > 1. This has 
been verified by comparing the 5^(24 fim)/S u (i band) 

2 see http:/ /irsa.ipac. caltech.edu/data/GOODS/docs/goods_dr3. html 



A Generalized Power-Law Diagnostic for Infrared Galaxies at z > 1 



3 



1000 



o 

O 

J, 

B 

CO 

Q 

x 

=5 



100 



1 1 1 1 1 1 1 1 1 

n 


i 1 i i i i 1 i i i i 1 i i i i 1 i i i i 


- 


o " 
- 


- ° n 


o ° 


r O 




\ 
\ 


o 


% °£ 






9 ULIRGs ° 




3 Oo ° 




o o 

§S e° ° 


Ipt 


s ^ ° o 








sPoN^ : 


LIRGs 


\ 

i 1 i i i i l\ i i i 1 i i i i 1 i i i i 



1 1.5 2 2.5 3 3.5 4 
Spectroscopic Redshift 

Fig. 1. — Flux densities S^(24/xm) vs. spectroscopic redshifts 
£spec for the 174 galaxies in the 24 /im sample with z S pec > 1- 
The dashed line, which divides the LIRG and ULIRG regimes, 
has been computed by considering the Bavouzet et al. (2008) 
monochromatic-to-total IR luminosity conversion, re-calibrated for 
galaxies with vL%^ > 10 10 L (see I Caput i~ et al. 2007). The k- 
correction factors are ba sed on a mixture of IR star-forming galaxy 
templates (Lagache et al. 2004) and empirical IR spectra. 



colour distributions of the current spectroscopic sample 
and the overall GOODS-S 24 /im sample with z > 1, 
with spectroscopic and photometric redshifts (Caputi et 
al. 2006a,b; 2007). 
The S u (24jam) flux densities versus spectroscopic red- 



shifts 



'spec 



of the 174 galaxies are shown in Fig. [TJ The 



dashed curve in this plot separates the LIRG and ULIRG 
regimes at bolometric IR luminosities Ljf> = 1O 12 L . 
To compute total IR luminosities, I have considered 
the Bavouzet et al. (2008) monochromatic-to-total IR 
luminosity conversion, re-calibrated for g alaxies with 
uL s^m > io^L© fsee ICaputi et al. 20071) . As in this 
work, the k-correction factors are b ased on a mixture o f 
IR star-forming galaxy templates (|Lagache et al. 2 004) 
and empirical IR spectra. Although, in a strict sense, 
this LIRG/ULIRG separation curve corresponds to star- 
forming galaxies, it is suitable for the purpose of illus- 
trating the appoximate regions of the LIRG and ULIRG 
regimes in the S V (2A jim)-z diagram. 

Independently, I made use of the publicly released 
GOODS-S IRAC maps a nd the SEXTRACTOR software 
(|Bertin fc Arnouts 19961 ) to extract a catalogue of IRAC 
sources at 3.6 and 4.5 /im. I measured the photometry at 
5.8 and 8.0 /im using SEXTRACTOR in dual-map mode. 
In order to construct a multi- wavelength optical/near- 
IR catalogue for the IRAC-selected sources, I searched 
for counterparts of these sources in the GOODS-S Very 
Large Telescope (VLT) ISAAC near-IR images, and the 
Hubble Space Telescope Advanced Camera for Surveys 
(ACS) maps, within a matching radius of 0.5 arcsec. 
The finally compiled catalogue includes photometry in 
11 bands: B, V, z, z, J, H, K 8 , [3.6], [4.5], [5.8] and [8.0]. 
All the magnitudes in this catalogue are total and have 
been corrected for galactic extinction. To obtain total 
magnitudes, I considered aperture magnitudes in circles 
of 4-arcsec (IRAC) and 2-arcsec (ISAAC and ACS) diam- 




Observed Wavelength (microns) 



Fig. 2. — Illustration of the technique applied here to identify 
galaxies with a significant SED power-law component. A power- 
law (dashed line) is subtracted from the original photometry of a 
galaxy (open circles), producing a new photometric dataset (filled 
circles). Stellar templates are then fitted to the original and power- 
law subtracted photometry (thin black and thick red lines, respec- 
tively). For each galaxy, 31 power-law subtractions with different 
combinations of a = 1.3, 2.0 and 3.0, and b = 0.0, 0.1, 0.2, . . . , 1.0, 
have been tested. In the example shown here, it can be seen that 
a power-law subtraction improves the overall SED fitting, partic- 
ularly allowing for a better agreement with the datapoints in the 
IRAC bands. 

eter, and applied the corresponding aperture corrections. 

Within the catalogue of IRAC-selected sources with 
multi- wavelength photometry, I identified all the 174 
24 /im galaxies with spectroscopic redshifts z spec > 1. 
The SED fitting analysis explained in next section is 
based on this photometry. A total of 22 galaxies (12.5% 
of the sample) are out of the field of view in all the ACS 
bands, or at least one of the ISAAC bands. For this mi- 
nority of sources, the SED fitting is based on 7,8,9 or 
10 bands, depending on the case. For all the remaining 
galaxies, the SED fitting is based on the 11-band pho- 
tometry. In case of a non-detection, a 2a upper limit has 
been considered for the flux density in the corresponding 
band. 

3. SPECTRAL ENERGY DISTRIBUTION ANALYSIS 

3.1. The Incidence of a Power- Law Component in the 
SEDs of Infrared Galaxies 

I considered that the flux density S v (A) of each galaxy 
at observed wavelengths A < 8 /im can be decomposed 
as: 



Si/ (A) = Sf elL (A) + 6^(8 /im) x 



A 



8 /im 



(i) 



where the two terms on the right-hand side correspond 
to a stellar component and a power-law component with 
spectral index a, respectively. The power-law term is 
normalised to the 8 /im flux density, so the constant b 
can vary between and 1. This decomposition is similar 
to that proposed by Hainline et al. (2011) to analyse the 
rest optical/near-IR SEDs of sub- millimetre galaxies at 

For each of the 174 galaxies in my 24 /im sample, 
I produced a suite of 31 photometric sets 5^ tell '(A) in 
the BVizJHK s [3.6] [4.5] [5.8] [8.0] bands, corresponding 
to different power-law subtractions. I considered that a 
could take three, and b ten possible values: a = 1.3,2.0 
and 3.0, and b = 0.0, 0.1, 0.2, 1.0. The case with 
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Fig. 3. — Minimum % 2 value obtained in the SED fitting of 24 /im 
galaxies with z S pec > 1 using stellar templates, after subtracting 
a power-law component to the photometry, vs. the original value 
obtained with no power-law subtraction. The dashed lines delimit 
the regions with A% 2 > 4.0 and 9.0, i. e. where the SED fitting 
with no previous power-law subtraction can be rejected with 2a 
and 3cr confidence, respectively. 



b = 0.0 and any a value corresponds to the original pho- 
tometry S^ telL (A) = Si/ (A). The three adopted a values 
are representative of the IR power-law indices that char- 
acterise different types of AGN, with different hot-dust 
components (e.g. Alonso-Herrero et al. 2003, 2006; Pol- 
letta et al. 2006; Honig et al. 2010). 

I performed the SED modelling on the 31 pho- 
tometric sets 5^ telL (A) of each galaxy, using a cus- 
tomised x 2 - m ini m isation SED-fitting code that incorpo- 
rates the 2007 version of the Bruzual & Chariot syn- 
thetic stellar template l i brary, with so lar metallicity 
(jBruzual fc Chariot 200l iBruzual 2QQ7h . These tem- 
plates correspond to different star formation histories: 
a single stellar population, and different exponentially- 
declining star- formation histories with r = 0.1 through 
5 Gyr. In all cases, the redshift of the galaxy has been 
fixed to the known z spec value. To account for inter- 
nal extinction, I convolved the stellar templates with 
the Calzetti et al. (2000) reddening law, allowing for 
0.0 < Ay < 3.0 with a step of 0.1. Stellar masses are 
obtained in the output of the same SED-fitting code. 

The technique is illustrated in Figure [2j once a power- 
law component is subtracted from the original photome- 
try (open circles), a new photometric dataset is obtained 
(filled circles). For each galaxy, I have tested 31 possible 
power-law subtractions (including the original photome- 
try), as explained above. I have then fitted stellar tem- 
plates to all these photometric sets, and compared the 
resulting minimum x 2 values. 

Figure [3] shows the minimum reduced x 2 value after 
power-law subtraction (i.e. the absolute minimum ob- 
tained considering all the different (a, b) combinations), 
versus the minimum reduced x 2 value obtained with the 
original photometry, for each galaxy. The points lying 
on the identity line correspond to galaxies for which the 
original photometry, without a power-law subtraction, 
produces the absolutely best fitting. The points below 
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Fig. 4. Same as Fig.ffl with the PLCGs highlighted. 



the identity line indicate the galaxies for which the SED 
fitting with stellar templates improves after the subtrac- 
tion of a power-law component. This is the case for the 
majority of the 24 /am galaxies. In particular, in 60 out 
of 174 cases (35%), the best SED-fitting solution on the 
original photometry, without a power-law subtraction, 
can be rejected with > 2a confidence. These 60 galaxies 
are, according to my definition, the PLCGs in the 24 /im 
sample with z spec > 1. 

Interestingly, around two-thirds of the PLCGs (42 out 
of 60) lie at z spec < 1.5 (Fig.HJ). This implies that - 47% 
of all the LIRGs and ULIRGs at 1.0 < z < 1.5 are 
PLCGs, i.e. they have a significant power-law compo- 
nent in their rest-frame near-IR SEDs. The percentage 
of PLCGs at 1.0 < £ sp ec < 1.5 is much higher than the 
known percentage of AGN among LIRGs and ULIRGs 
at these redshift s, and raises the question of whether the 
power-law component is actually related to an AGN pres- 
ence in all these cases. I investigate this issue further in 
Section At z > 1.5, instead, the fraction of PLCGs 
among 24 /am galaxies (the vast majority of which are 
ULIRGs) is of only - 21%. 

I repeated the same SED modelling using the 2003 ver- 
sion of the Bruzual & Chariot stellar template library, 
and found very similar results: ~ 38% of the 24 /xm 
sources with z spec > 1 are classified as PLCG. The result- 
ing minimum x 2 values are very similar to those obtained 
with the 2007 templates (the median of the minimum x 2 
differences is 0.009). All the following analysis will be 
based on the 2007 template library run, but no conclu- 
sion in this paper would change if I adopted the results 
obtained with the 2003 templates. 

The host galaxy properties, as determined from the 
best-fitting stellar templates, are similar for PLCGs and 
non-PLCGs. For both sub-samples, the best-fitting star 
formation histories correspond to single stellar popula- 
tions or r models with r < 0.1 Gyr in > 80% of cases. 
After the power-law subtraction, these models provide 
the best fitting for 93% of the PLCGs. The best-fitting 
extinction values for PLCGs (before power-law subtrac- 
tion) and non-PLCGs are also very similar: the median 
and r.m.s. of the distributions are Ay = 1.4 ± 0.6 and 
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Fig. 5. — Best-fitting extinction Ay vs. spectroscopic redshifts 
for the 60 PLCGs, classified according to their spectral index a. 
The Ay values obtained before and after power-law subtraction are 
shown in each case (small open circles and larger coloured circles, 
respectively). The Ay — z spe c values of the non-PLCGs are shown 
as a reference in all panels (dots). 



Ay = 1.4 ± 0.7, respectively. After power-law subtrac- 
tion, the PLCGs have Ay = 1.0 ± 0.6. 

Figure [5] shows the extinction Ay versus spectroscopic 
redshifts for the 60 PLCGs, classified according to their 
power-law indices a. The non-PLCGs are also shown as 
a reference. A bit more than a half of the PLCGs (55%) 
are characterised by a power-law component with a = 
3, while the remaining PLCGs have lower best-fitting 
power-law indices, i.e. a = 1.3 or 2. These indices are 
related to the temperature of the hot-dust component in 
the galaxy, which is higher for lower a values (cf. Section 

The effect of a decrease in the best-fitting Ay value 
is especially evident for the a = 1.3 PLCGs: a shallow 
power-law component can mimic the effect of additional 
reddening. After power-law subtraction, the majority of 
the a — 1.3 and 2 PLCGs have Ay < 1, while most of 
the a = 3 PLCGs have Ay > 1. 

The normalisation factor of the power-law component 
is b > 0.50 for virtually all the PLCGs, which means that 
at least 50% of the 8 /im flux density is in the power- 
law component. For around a half of the PLCGs, this 
contribution is > 80%. This confirms that the PLCG 
definition adopted here, based on the > 2a improvement 
of the SED fitting, truly selects sources for which the 
power-law component makes an important contribution 
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Fig. 6. — Same plot as in Fig. [3] but for the control sample of 
247 IRAC galaxies with z spec > 1 which are not 24 /xm-detected. 

to the rest-frame near-IR light of the galaxy. 

For galaxies at z ~ 1.4, the emission from polycyclic 
aromatic hydrocarbons (PAHs) at rest A rest ~ 3.3 /im 
will enter the observed 8 jam band. Although this is a 
relatively faint PAH emission feature, it can potentially 
mimic the effect of a steep power-law component in the 
galaxy SED. To assess the importance of this effect on 
the PLCG selection, I have performed the SED fitting on 
additional photometric sets for each galaxy, with power- 
law subtractions corresponding to larger spectral indices, 
namely, a = 3.5 and 4. Within the total sample, only 
13 galaxies appear to be PLCGs with such large best- 
fitting spectral index. However, their redshifts are not 
concentrated around z ~ 1.4 or any other specific red- 
shift, so the steeper a values are unlikely the effect of 
emission features. In fact, 11 out of these 13 galaxies 
have been recognised as PLCGs with a = 3 in my orig- 
inal analysis. Only two galaxies appear as new PLCGs: 
one at z spec = 2.810 (corresponding to an X-ray luminous 
AGN), and another one at z spec = 1.550 (for which the 
3.3 /im PAH emission could be partly within the 8 /im- 
filter wavelength coverage). These results indicate that: 
a) the 3.3 /im PAH emission plays a minor role in the 
overall IR SED for most IR galaxies; b) considering a > 3 
values for the power-law subtraction produces a virtually 
negligible effect on the identified PLCG sample. 

To investigate whether the high fraction of PLCGs 
is truly a characteristic of the 24 /im galaxy sample, I 
performed a similar SED fitting analysis on a control 
sample of 247 galaxies selected from the IRAC cata- 
logue described in Section ^2j also with secure £ spe c > 1, 
but which are not 24 /im-detected (i.e. S v {2Ajim) < 
80/iJy). These galaxies have been selected to have a 
similar stellar mass distribution as the 24 /im galaxies 
(Msteii. > 10 10 M o , based on the SED modelling of the 
original photometry). For each of these galaxies, I have 
produced 31 sets of photometric variations in the same 
manner explained above, and fitted the resulting SEDs 
with stellar templates in all cases. The comparison of 
minimum y 2 values with and without power-law sub- 
traction is shown in Fig. [6] 
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Fig. 7. — Distribution of corrected-to-original stellar mass ratios for the PLCGs, classified according to their best-fitting power-law index 
a. The corrected stellar masses are those derived from the SED modelling after the best-fitting power-law subtraction. The vertical, dashed 
lines indicate the median correction factors derived from each distribution. 

For the control sample, the percentage of galaxies that 
significantly improve their SED fitting with stellar tem- 
plates after a power-law subtraction is a factor of two 
smaller than for the 24 /xm galaxies: a pure stellar SED 
fitting can be discarded with > 2a confidence in only 
18% of cases (44 out of 247 galaxies). These 44 PLCGs 
in the control sample have best-fitting a values shared 
in similar proportions as for the 24 /im-detected PLCGs: 
50% have best-fitting a = 3, while the other 50% have 
a = 1.3 or 2. 

These results obtained on the control sample indi- 
cate that the incidence of PLCGs among mid-IR-selected 
galaxies is twice more important than among other sim- 
ilarly massive galaxies at z > 1. 



3.2. The Effect of a Power-Law Subtraction on the 
Derived Stellar Masses 

Figure [3 shows the distribution of corrected-to-original 
stellar mass ratios for the 60 24 /xm-selected PLCGs, 
classified according to their best-fitting power-law in- 
dex. The corrected stellar masses are those derived from 
the SED modelling with stellar templates after the best- 
fitting power-law subtraction. 

The inspection of these distributions shows that the 
correction on the stellar masses is non-trivial: the exact 
correction needs to be derived on a case-by-case basis. 
However, for >80% (90%) of the PLCGs the correction 
is within a factor of ~ 1.5 (2.0). The widest correction- 
factor distributions are those corresponding to power-law 
spectral indices a = 1.3 and 2, for which the stellar mass 
ratios range between < 0.1 and ~ 1.8. The median 
values of these stellar mass ratios are 0.67 and 0.85, re- 
spectively. 

The corrections on the stellar masses are the smallest 
for the case a = 3 (the median of the ratio distribution 
is 0.92). This behaviour is expected, as the subtraction 
of such a steep power-law has a rapidly declining effect 
on the photometry at wavelengths shorter than 8 /mi, 
i.e the overall SED shape is not substantially modified 
except for a single photometric point. 

It may a priori seem surprising that the SED fitting 
for some galaxies yields larger stellar masses after the 
PL subtraction. This effect is produced because the 
best-fitting stellar template properties (star formation 
history, age and extinction) are usually different for the 
PL-subtracted photometry to those for the original pho- 



tometry, and, in some cases, the resulting overall nor- 
malisation factors (i.e. mass-to-light ratios) are larger. 

For 5 out of the 60 PLCGs, all with best-fitting 
a = 1.3, the corrected stellar mass is less than 1/10 of 
the stellar mass derived from the original photometry. 
These cases correspond to "pure power-law sources" , and 
all except one of them are identified as X-ray luminous 
AGN with L x > 10 43 erg s" 1 in the Chandra 4 Ms cata- 
logue for the GOODS-S field (jXue et al. 201 If ). The only 
source that is not X-ray detected is the highest-redshift 
of these pure power-law sources, with S v (2Ajim) « 
286 /iJy and z spec = 3.356, so it is likely an intrin- 
sically luminous AGN as well. For these sources, a 
stellar mass calculation based on the original photom- 
etry should clearly be avoided. But, as said above, 
in all other cases the correction to the stellar mass is 
within a factor of two. This result validates the ap- 
proach of obtaining approximate stellar mass estimates 
for e.g. low-luminosity AGN, using the original un- 
corr ected photometry, that is followed by som e authors 
(e.g. I Silverman et al. 20091 : iRosario et al. 20l3[ ). 

4. THE NATURE OF THE PLCGS 

4.1. Independent AGN Diagnostics: X-ray Detections 

In order to investigate whether the PLCGs actually 
host AGN, I searched for counterparts of these sources 
in the publicly available, ultra-de ep Chandra 4Ms X-ray 
catalogue for the GOODS-S field (jXue et al. 201 ID . This 
diagnostic is very useful to assess to which extent the 
PLCG criterion can be used to identify AGN in the gen- 
eral case that such deep X-ray catalogues are not avail- 
able (to date, X-ray maps of such a depth have only been 
obtained for the GOODS-S field). 

Fig. [8] shows the S V (2A jim)-z diagram for the 24 jam 
galaxies in the z spec > 1 sample that are identified 
with X-ray sources in the Chandra 4Ms catalogue. The 
PLCGs are highlighted with couloured circles. There 
are a total of 59 24 /im galaxies identified with an X-ray 
source, 48 out of which correspond to X-ray-classified 
AGN, and 11 are X-ray 'normal galaxies'. The X- 
ray classification is based on the intrinsic X-ray lumi- 
nosities, X-ray colours (hardness rat ios), and derived 
neutral gas co lumn densities (e.g. Bau er et al. 20041 : 
Kue et al. 201 ID . Although the X-ray-based AGN classi- 
fication is among the most secure, the criteria typically 
adopted to classify X-ray sources could still result in some 
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Fig. 8. — Flux densities 5^(24/xm) vs. spectroscopic redshifts 
Zspec for the 24 jum galaxies at z S pec > 1 that are identified with 
X-ray sources in the Chandra 4Ms catalogue. Asterisks and crosses 
correspond, respectively, to AGN and normal galaxies, classified 
according to their X-ray properties. All PLCGs are shown with 
couloured circles. 

low-luminosity or highly-obscured AGN being classified 
as normal galaxies (see Xue et al. 2011 and references 
therein). 

The analysis of Fig. [8] shows the following. 
• At z < 1.5, 77% of the X-ray AGN are recognised as 
PLCGs, which indicates that the PLCG criterion pro- 
duces a high completeness in the identification of AGN 
at these redshifts. And, interestingly, among the X-ray- 
classified 'normal galaxies', the ratio of PLCGs is very 
similar (~ 78%). So, either these galaxies also host an 
AGN, or there is another mechanism that gives rise to the 
X-ray emission, and is also responsible for the hot-dust 
component manifested in the PLCG classification. 

In starburst galaxies, supernova explosions can pro- 
duce hot-gas bubbles with significant X-ray emission. 
The resulting shock waves can heat the surrounding 
dust grains up to sublimation temperatures T > 1000 K. 
However, if the X-ray luminosities of X-ray galaxies were 
only a product of supernova explosions, they would corre- 
late with the IR lumi nosities, as both trace the on-going 
star formation rates ([Mas- Hesse et al. 2008[ ). This is not 
the case for the PLCGs at z < 1.5, for which two sources 
can differ in more than a factor of three in IR luminos- 
ity, but both be classified as X-ray normal galaxies with 
similar X-ray luminosities. This becomes clear by in- 
spection of Fig. 2J where some bright 24 /im sources at 
1.0 < z < 1.5 are neither PLCGs nor X-ray galaxies at 
all. Thus, the SED power-law excess of these PLCGs is 
not simply tracing high star formation rates. 

At z < 1.5, it is interesting to note the different distri- 
butions of best-fitting spectral indices a among the X-ray 
AGN and normal galaxies that are recognised as PLCGs: 
9 out of 13 (69%) X-ray AGN / PLCGs have a = 1.3 or 2, 
while only 2 out of 7 (29%) of the X-ray galaxy / PLCGs 
have a = 1.3 or 2. Or, equivalent ly, only 31% of the X- 
ray AGN classified as PLCGs have a = 3, while 71% of 
the X-ray galaxies classified as PLCGs have a = 3. I 
will carry on this discussion in Section ^4.3[ where I will 



argue that the spectral index a is a good discriminator 
of PLCGs of different nature. 

• Only 42% of the X-ray-detected AGN at z > 1.5 
are classified as PLCGs. This effect could be due to 
a real difference in the properties of X-ray AGN at dif- 
ferent redshifts: e.g. at z > 1.5 the light of the host 
galaxies may dominate the rest near-IR SED, making a 
power-law component to appear less significant. Alter- 
natively, the lower percentage of AGN incidence could 
be the consequence of a subtle k-correction effect, as 
the observed 8 /im photometry only traces rest wave- 
lengths A rest < 3/im at z > 1.5, while it traces up to 
rest A r est ~ 4 /im at z = 1. 

To investigate these two possible explanations, I con- 
sider those 24 /im galaxies at z ~ 2 that form part of the 
Fadda et al. (2010) IR galaxy sample. Pozzi et al. (2012) 
have studied the mid-/far-IR SED decomposition of these 
galaxies, and found that nine out of 24 of them have 
a significant AGN component. Only three out of these 
nine sources are recognised as PLCGs within my sample. 
This result suggests that data beyond observed 8 /im are 
necessary to reveal more effectively an AGN presence in 
galaxies at z ~ 2. 

In Section £14. 4[ I will show that very hot dust tem- 
peratures (Tdust 1300 K) are necessary to generate 
a power-law SED component detectable in the IRAC 
bands at z > 1.5. Su ch hot dust is typically found i n 
AGN inner dusty tori (jBarvainis 1987MMor et al. 2012( 1. 
but some AGN are characterised by lower dust tem- 
peratures, or the inner torus may be less exposed (e.g. 
iSchartmann et al. 2005[ ). This will restrict the fraction 
of AGN that can be identified with my PLCG criterion, 
which is based on SED fitting up to observed 8 /im, at 
z~2. 

• Around a half of the PLCGs at 1.0 < z < 1.5 are not 
X-ray detected. The mean optical/near-IR SED prop- 
erties of these PLCGs are not significantly different to 
those of X-ray-detected PLCGs classified as X-ray galax- 
ies: e.g. the median extinction is Ay = 1.0 for the X-ray 
undetected, against Ay = 0.8 for the X-ray detected. 
Their IRAC colours are also similar (see Section $5}. 
The PLCG classification, and the X-ray detection, do 
not follow a simple correlation with the host galaxy IR 
luminosities: there are X-ray detected PLCGs both in 
the LIRG and ULIRG regimes. But all the X-ray unde- 
tected PLCGs are LIRGs, which suggests that the X-ray 
non-detection is due to the sensitivity limit of the X-ray 
surveys. 

4.2. PLCGs Non-Detected in X-Rays: Stacking Analysis 

To probe whether the X-ray non-detection of half of the 
PLCGs at 1.0 < z < 1.5 is simply due to the sensitivity 
limit of the Chandra 4Ms X-ray survey, I performed a 
stacking analysis of these sources in X rays. 

Fig. [9] shows the stacked X-ray images of the PLCGs 
that are not individually detected in X rays at 1.0 < z < 
1.5. In the hard X-ray band (2.0-8.0 keV), the stacking 
does not produce any detection, while in the soft band 
(0.5-2.0 keV) a detection with a > 6<r significance is ob- 
tained. This stacked signal confirms that these sources 
are X-ray emitters, but their fluxes are below the sensi- 
tivity limits of the 4Ms Chandra survey. 

Nuclear activity cannot be discarded as an agent partly 
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Fig. 9. — Stacked X-ray images of the 1.0 < z sp ec < 1.5 PLCGs 
that are not individually detected in X rays. Left: 0.5-2.0 keV; 
right: 2.0-8.0 keV. Each stamp covers 11 x 11 arcsec 2 . 

responsible for the soft X-ray emission in these sources, 
but it is probably not the main cause: on average they 
have log-m(/x ///?.) ~ —2.7, too low a value for an AGN 



classification (jXue et al. 20111 ). 

Around two thirds of the stacked PLCGs have spectral 
index a = 3 which, as I discuss below, should mainly 
correspond to star-forming galaxies with hot interstel- 
lar media (ISM) at 1.0 < z < 1.5. The stacking of 
the a = 3 PLCGs alone yields a 4.8<r detection in the 
(0.5-2.0 keV) band, which corresponds to an average 
flux density fx ~ 1.9 x 10 -1T ergcm -2 s _1 (after cor- 
rection for galactic extinction; see Stark et al. 1992). At 
the median redshift of these sources z me( ± = 1.13, this 
corresponds to a luminosity Lx ~ 1.3 x 10 41 ergs _1 , 
assuming a photon index T = 2, which is appro priate 
for star- forming galaxies (e.g. iLehmer et al. 20 08). Us- 
ing the relation calibrated by Mineo et al. (2013), I find 
that this corresponds to an instantaneous star formation 
rate SFR w (37 ± 4) M yr _1 . This is significantly lower 
than the SFR derived from the average total IR lumi- 
nosities of these sources, which is SFR « 9OM yr _1 , 
as obtained using the SFR-Lir relation given by Ken- 
nicutt (1998). This difference indicates the presence of 
internal absorption. To reconcile both SFR values, I find 
that these galaxies should have an average column den- 
sity N H ~ 4.5 x 10 21 cm -2 . 

In any case, as commented above, one has to be careful 
by considering the PLCG X-ray luminosity as a tracer of 
the instantaneous star formation rate, as many of these 
galaxies are likely composite star- forming/ AGN systems. 
Also, the X-ray emission may be the result of energetic 
processes in the galaxy ISM that are not simply related 
to the on-going star formation rate. Some of these effects 
probably account for the dispersion observed in the Lx~ 
SFR relation. 

4.3. Morphologies of PLCGs and X-Ray Galaxies at 
l.0<z< 1.5 

Further clues about the nature of PLCGs, particu- 
larly at 1.0 < z < 1.5, can be obtained by comparing 
their morphologies. Figure fTUI shows HST/ ACS z-band 
postage stamps of 15 PLCGs at 1.0 < z < 1.5, grouped in 
three different categories: X-ray detected sources, classi- 
fied as AGN and Galaxies, and PLCGs that are not in- 
dividually X-ray detected. The three classes display ob- 
jects with different morphologies, including galaxies with 
compact and regular matter concentrations, and galaxies 
with very disturbed morphologies. The latter are seem- 
ingly the product of galaxy mergers, or are characterised 
by other violent processes, such as the radial ejection of 
material (see the X-ray galaxy at the right end of the 



middle row). 

The striking conclusion that can be extracted from 
FigJlO] is the trend between the PLCG morphology and 
the value of the best-fitting spectral index a: those 
PLCGs with a = 1.3 or 2.0 show quite a regular mor- 
phology, while the PLCGs with a = 3.0 are irregular and 
show signs of different kinds of galaxy interactions. 

This trend is independent of the X-ray detection and 
classification. This fact suggests that the spectral index 
a contains information about the nature of the source: 
in the PLCGs with a = 1.3 or 2.0, the SED power-law 
component is very likely the signature of nuclear activ- 
ity, while in the a = 3 PLCGs the power law appears 
to be the consequence of a hot ISM produced by the vi- 
olent galaxy dynamics. Shocks of hot gas can compress 
and heat up the ISM dust to several hundred Kelvin, 
and even up to sublimation temperatures (T > 1000 K; 
lOsterbrock fc Ferland 2006( 1. This hot dust, not related 
to an AGN origin, will also manifest itself as a power-law 
component in the rest near-IR galaxy SED. 

Of course, in some cases, an AGN and a host galaxy 
with a disturbed ISM may co-exist, producing X-ray con- 
firmed AGN that are PLCGs with a = 3 (upper right 
stamps in Fig. [TUJ see also Kocevski et al. 2012 for a 
discussion on AGN morphology). 

4.4. Estimated Temperatures of the PLCG Hot-Dust 
SED Components 

I estimated the dust temperatures necessary to pro- 
duce a power-law component in the galaxy SEDs, as- 
suming a greybody function for the dust emission: 
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where the second factor on the right-hand side of equa- 
tion ([2]) corresponds to the black-body emission at tem- 
perature T, and h and k are the Planck's and Boltz- 
mann's constants, respectively. I assumed /3 = 1.5, as it 
is common in the literature. 

Fig. [11] shows the observed flux density ratios 
S' Z y(8/im)/S Z y(4.5 jam) produced by greybody emitters at 
different temperatures, as a function of redshift. The hor- 
izontal solid lines indicate the fixed, observed flux density 
ratios produced by power-laws with the three spectral in- 
dices considered in this paper, i.e. a = 1.3, 2.0 and 3.0. 

Two main results can be derived from this simple dia- 
gram. Firstly, to generate a power-law component with 
a = 1.3 or 2 at 1.0 < z < 1.5, temperatures as high as 
900 — 1000 K are needed. Instead, for a power-law with 
spectral index a = 3.0, temperatures T ~ 700 — 800 K 
are sufficient. This may explain why most of the X-ray 
AGN, for which dust temperatures close to the subli- 
mation point are expected, are associated with PLCGs 
with spectral indices a = 1.3 or 2, at 1.0 < z < 1.5. The 
PLCGs with spectral index a = 3.0 can be explained 
with a shock- heated ISM, with no need of dust getting 
close to the sublimation point (the exact dust sublima- 
tion temperature depends on the dust composition, but 
it is typically T > 1000 K). The existence of hot inter- 
stellar dust in IR galaxies has been previously discussed 
in the literature. Particularly, Lu et al. (2003) found that 
a colour dust temperature of ~ 750 K can account for the 
near-IR excess observed in some local IR normal galaxies, 
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Fig. 10. — HST/ ACS z-band postage stamps of 15 PLCGs at 1.0 < z < 1.5. Each stamp shows an area of 5 x 5arcsec 2 centred at the 
source. The best-fitting spectral index a is indicated in each case. The solid lines separate the cases with a = 3.0, which have a much more 
disturbed moroholoerv than the PLCGs with smaller soectral indices a = 1.3 and 2.0. 

to be more relia ble for AGN identif ication among fainter 
IR sources (e.g. iDonlev et al. 20121 ). Here I analyse how 
the PLCGs, and other 24 /mi galaxies at z spec > 1, are 
classified in an IRAC colour-colour diagram. 

Fig. [12] shows the IRAC colour-colour diagram for all 
the 24 /im galaxies at z spec > 1, with the PLCGs and 
X-ray sources highlighted. The colours shown here are 
overall colours, i.e. there is no separation of stellar 
and power-law SED components. One can see that the 
range of overall IRAC colours displayed is quite wide, 
with some 24 /im galaxies having much redder IRAC 
colours than others. PLCGs, in particular, display a 
wider range of colours than the non-PLCGs, as there 
are very few non-PLCGs with very red IRAC colours 
log 10 (S 8 .o/S 4 .5) > 0.2 and log^Ss.s/Sa.e) > 0.1. For 
PLCGs, I find no correlation between the IRAC colours 
and the spectral index of the SED power-law component. 

The solid lines in Fig. [T2l delimit the revised colour re- 
gion for AGN selection proposed by Donley et al. (2012). 
More than 80% of the 24 /im galaxies that lie within that 
region are PLCGs or X-ray AGN, indicating that indeed 
the proposed colour selection is highly reliable. However, 
these colour criteria are too strict and miss a substan- 
tial fraction of PLCGs, many of which are X-ray con- 
firmed AGN. The incompleteness of this colour selection 
was discussed by Donley et al. (2012), who argued that 
their criteria quickly lose completeness for AGN with X- 
ray luminosities Lx < 10 44 ergs _1 . Within the current 
sample, all the X-ray-detected sources in the Donley et 
al. (2012) region have L x > 5 x 10 43 ergs -1 . 

There are three galaxies that lie within the Donley et 
al. AGN colour-colour region, but which are neither iden- 
tified as PLCGs, nor X-ray detected. One of them is a 
^spec = 1.045 galaxy for which the SED fitting with a 
power-law subtraction just lies below the 2a criterion 
for PLCG identification. So, this galaxy could still have 
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Fig. 11. — Observed flux density ratios 5^(8.0 /xm)/Si/(4.5 /xm) 
for greybody emission at different temperatures, as a function of 
redshift (dashed lines). The horizontal, solid lines indicate the 
fixed, observed ratios produced by power-law components with dif- 
ferent spectral indices a. 

similarly to what I find here for LIRGs at 1.0 < z < 1.5. 

Secondly, at redshifts z > 1.5 — 2.0, producing a power- 
law component in the observed IRAC SED requires hot 
dust with temperatures T > 1300 K, for any of the spec- 
tral indices. This mostly explains why the incidence of 
PLCGs is lower at z > 1.5, as only sources with such an 
extremely hot dust component will be classified as PLCG 
with the criterion adopted in this paper. 

5. IRAC COLOUR-COLOUR DIAGNOSTIC 

The use of Spitzer/IRAC colour-colour diagrams has 
proven to be quite s uccessful for identifying AGN amon g 
IR-selected galaxies ([Lacv et al. 20041 : IStern et al.~2 005). 
Over the years, these colour selections have been refined 
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a buried AGN. The other two galaxies are at redshifts 
^spec = 2.442 and 2.577. For these, the stellar bump cen- 
tred at rest-frame ~ 1.6 /am is within the observed 5.8 /am 
filter, producing a red S^.Sfim/Ss.Qfim colour. These two 
galaxies are close to the ^.o^m/^.s^m lower boundary 
of the Donley et al. colour wedge, so their presence 
within the wedge may simply be due to photometric scat- 
tering (note the size of the error bars). 

Galaxy colour selections are attractive for their sim- 
plicity, but no criterion has been found to safely identify 
the bulk of the AGN population at any redshift. The 
PLCG criterion proposed here offers a suitable alterna- 
tive, which is particularly successful at identifying AGN 
of different luminosities at 1.0 < z < 1.5. Note that it 
would not be enough to extend the Donley et al. (2012) 
to bluer colours, because most of the PLCGs display sim- 
ilar IRAC colours as other IR sources which do not have 
a significant power-law component in their SEDs, so an 
extended colour selection would produce a highly con- 
taminated AGN sample. 

6. SUMMARY AND CONCLUSIONS 

In this paper I have studied the rest-frame 
UV/optical/near-IR (A rest < 3/im) SEDs of 174 24 ^in- 
selected galaxies with secure spectroscopic redshifts 
z spec > 1, analysing explicitly the presence and signif- 
icance of a hot-dust, power- law component. I modelled 
the residual light produced after subtracting different 
possible power-law components, using stellar templates, 
on 11 broad-bands from the 5-band through 8 jam. The 
subtracted power-law components are characterised by 
three different spectral indices a = 1.3, 2.0 and 3.0, and 
different normalisation weights. 

I found that in 60 out of 174 (35%) cases, the SED fit- 
ting with stellar templates without any previous power- 
law subtraction can be rejected with > 2<j confidence, i.e. 
35% of the 24 jam galaxies at z > 1 are characterised by 
a significant power-law component in their rest near-IR 
SED. I referred to these galaxies as PLCGs. This high 
percentage of PLCGs at z > 1 appears to be inherent 
to the IR bright galaxy population. A similar analy- 
sis performed on a control sample of 247 IRAC-selected 
galaxies with z spec > 1 and similar stellar masses as the 
24 /im galaxies, but which are not 24 /im-detected, shows 
that the incidence of PLCGs is of only 18%. 

The incidence of PLCGs among 24 jam galaxies is much 
higher at 1.0 < z < 1.5 (47%) than at z > 1.5 (21%). 
The lower incidence at z > 1.5 is explained by the 
fact that the presence of very hot dust (Td us t > 1300 K) 
is necessary to have a significant power-law component 
manifested at observed wavelengths A < 8 /im. There- 
fore, only sources with the most extreme dust conditions 
would be recognised with the PLCG criterion at these 
higher redshifts. 

The impact of a significant power-law component in the 
rest near-IR SED on the derived galaxy stellar masses 
needs to be analysed on a case-by-case basis, but it is 
small in general. The correction for stellar masses is 
within a factor of 1.5 in > 80% of cases. Only in a small 
percentage (< 10%) of the PLCGs, basically all corre- 
sponding to powerful AGN, the correction to the stellar 
mass is so large that a stellar mass calculation based on 
the original photometry should clearly be avoided. 
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IRAC colour-colour diagram for all the 24 sources 
spec > 1 (circles). As in previous figures, PLCGs are high- 
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lighted in colour. Asterisks and crosses indicate, respectively, X- 
ray-classified AGN and galaxies. The solid lines delimit the colour 
region for AGN selection proposed by Donley et al. (2012). Typical 
error bars on the colours are shown in the bottom right corner of 
the plot. 

A key issue is understanding whether the PLCGs are 
truly related to the presence of an AGN, particularly 
at 1.0 < z < 1.5, where the PLCG incidence is very 
high. To investigate this, I studied different properties 
of these sources. Firstly, I looked for identifications of 
the PLCGs in the ultra-deep 4Ms X-ray catalogue for 
the CDFS. I found that 77% of the X-ray AGN within 
my 24 jam sample with 1.0 < z spec < 1.5 are classified as 
PLCGs. This indicates that the PLCG criterion provides 
a high completeness in selecting AGN at these redshifts. 
Interestingly, 78% of the X-ray normal galaxies within 
my 24 /mi sample with 1.0 < z spec < 1.5 are also PLCGs. 
But their PLCG classification or X-ray luminosities do 
not correlate with the 24 /am luminosities, so the X ray 
emission or the power-law SED excess is not simply a 
tracer of high star formation rates. 

A stacking analysis of the PLCGs that are not indi- 
vidually X-ray detected at 1.0 < z < 1.5 reveals that 
they have on average similar properties as X-ray-detected 
normal galaxies. The stacked X-ray images yield a sig- 
nificant detection in the soft X-ray band (0.5 — 2.0 keV), 
and no signal in the hard X-ray band (2.0 — 8.0 keV). 
Around two thirds of the stacked PLCGs have a best- 
fitting power law with a = 3. 

The visual inspection of optical HST/ ACS images has 
proven very useful to shed light on the nature of the 
PLCGs at 1.0 < z < 1.5. I found a clear trend relat- 
ing the PLCG morphology and the spectral-index a of 
the best-fitting SED power-law component: PLCGs with 
a = 1.3 or 2.0 look regular and typically have a nuclear 
matter concentration. Instead, the PLCGs with a = 3.0 
typically show irregular morphologies, indicating a dis- 
turbed galaxy dynamics, which in some cases suggest 
the presence of galaxy mergers. This morphology trend 
is independent of the PLCG X-ray detection and classifi- 
cation. The a = 3.0 power-law component appears then 
to be a signature of quite extreme ISM conditions, where 
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the interstellar dust is heated to temperatures T > 700 K 
by the gas shocks produced in the disturbed ISM. 

Of course, this does not preclude the a = 3.0 PLCGs 
to also host AGN. Actually, the presence of an AGN 
is confirmed in some of them by the X-ray data. But 
an AGN presence does not appear to be necessary to 
produce a significant a = 3.0 power-law component. For 
the PLCGs with a = 1.3 or 2.0, instead, higher dust 
temperatures are necessary T > 900 K. This property, 
along with the regular morphologies, take to the direct 
conclusion that these sources must host an AGN. One 
could speculate that the a = 3.0 PLCGs may all be on 
the way to form an AGN as well, but this hypothesis 
cannot be tested with the current data. 

Considering that all the PLCGs with a = 1.3 or 2.0 in 
my sample contain AGN, and adding all the other X-ray 
classified AGN within my 24 /im sample, I obtain that a 
total of 30% of the 24 /xm sources at 1.0 < z < 1.5 host 
AGN. This should be considered as a lower limit, as some 
of the a = 3.0 PLCGs that are not individually X-ray de- 
tected could also contain AGN. If all of them had AGN, 
that would give an upper limit of ~ 52% on the AGN 
fraction among LIRGs and ULIRGs at 1.0 < z < 1.5. 
Very recently, different works have concluded on a higher 
fraction of AGN among IR galaxies than previously 
known at z < 1, obtaining percentages of 37- 50% (e.g. 
lAlonso-Herrero et al. 20121 : Uuneau et al. 20131 ). The re- 
sults of my study are in line with such conclusions. At 
z > 1.5, a similar criterion as that used at z < 1.5 would 
indicate that ~ 40% of the 24 /am galaxies contain an 
AGN, although this figure is largely based on the X-ray 
detections, rather than the PLCG incidence, which is rel- 
atively low (for a comparison, Pozzi et al. 2012 derived a 
fraction of 35% based on the mid-/far-IR SED analysis 
of 24 ULIRGs at z ~ 2). 

In this work, I have calibrated the PLCG analysis tech- 
nique making use of galaxies with secure spectroscopic 
redshifts. A complete analysis of the applicability of this 
technique with the simultaneous determination of pho- 
tometric redshifts (%>hot) will be presented in a future 
paper. A preliminary run of my SED fitting code leaving 
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the redshift as a free parameter indicates that, in fact, the 
consideration of a power-law subtraction produces some 
improvement in the overall £ p hot-^spec comparison, reduc- 
ing the number of catastrophic outliers. This suggests 
that the identification of PLCGs can be done, and the 
power-law subtraction is actually recommendable, when 
attempting a %>hot determination. 

Although less straightforward than a simple IRAC 
colour-colour selection, the PLCG identification criterion 
introduced here is much more complete to select AGN 
of different luminosities, especially in the redshift range 
1.0 < z < 1.5. Keeping only those PLCGs with spectral 
index a = 1.3 and 2 should result in the best compro- 
mise between completeness and reliability, as most X-ray 
normal galaxies and X-ray non-detections are among the 
PLCGs with a = 3. In the absence of very deep X-ray 
data, the PLCG selection with spectral index segrega- 
tion offers the most efficient method to identify the AGN 
presence in any large galaxy sample. 
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